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Abstract 

In  this  paper,  an  approach  for  fault  diagnosis  of  hybrid 
dynamic  systems  (HDS),  in  particular  discretely  controlled 
continuous  system,  is  proposed.  The  goal  is  to  construct  a 
decentralized  diagnosis  structure,  able  to  diagnose 
parametric  and  discrete  faults.  This  approach  considers  the 
system  as  composed  of  a  set  of  interacted  hybrid 
components  (HCs).  Each  HC  is  composed  of  a  discrete 
component  (Dc),  e.g.  on/off  switches,  with  the  continuous 
components  (Ccs),  e.g.  capacitors,  whose  continuous 
dynamic  behavior  is  influenced  by  the  Dc  discrete  states.  A 
local  hybrid  diagnosis  module,  called  diagnoser,  is 
associated  to  each  HC  in  order  to  diagnose  the  faults 
occurring  in  this  HC.  In  order  to  take  into  account  the 
interactions  between  the  different  HCs,  local  diagnosis 
decisions  are  merged  using  a  coordinator.  The  latter  issues  a 
final  decision  about  the  origin  of  the  fault  and  identifies  its 
parameters.  The  advantage  of  the  proposed  approach  is  that 
local  hybrid  diagnosers  as  well  as  the  coordinator  are  built 
using  local  models.  The  proposed  approach  is  applied  to 
achieve  the  decentralized  diagnosis  of  discrete  and 
parametric  faults  of  power  electronic  three-cell  converters. 

1.  Introduction 

1.1  Basic  definitions  and  motivation 

A  fault  can  be  defined  as  a  non-permitted  deviation  of  at 
least  one  characteristic  property  of  a  system  or  one  of  its 
components  from  its  normal  or  intended  behavior.  Fault 
diagnosis  is  the  operation  of  detecting  faults  and 
determining  possible  candidates  that  explain  their 
occurrence.  Most  of  real  systems  are  hybrid  dynamic 
systems  (HDS)  (Zaytoon,  2001),  (Arogeti  et  al .,  2010)  in 
which  the  discrete  and  continuous  dynamics  cohabit. 
Therefore,  fault  diagnosis  of  HDS  must  deal  with  the 
evolution  of  continuous  dynamics  in  each  discrete  mode  in 
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order  to  construct  a  diagnosis  module  (called  diagnoser) 
able  to  diagnose  parametric  and  discrete  faults.  Parametric 
faults  affect  the  system  continuous  dynamics  and  are 
characterized  by  abnormal  changes  in  some  system 
parameters;  whereas  discrete  faults  affect  the  system 
discrete  dynamics  and  are  considered  either  as  the 
occurrence  of  unobservable  events  and/or  reaching  discrete 
fault  modes.  In  both  cases,  they  entail  unpredicted, 
abnormal,  change  in  the  system  configuration.  Therefore 
faults  may  be  modelled  in  HDS  by  introducing  parameters 
into  the  system  model,  explicit  fault  events  or/and  fault 
modes. 

Discretely  controlled  continuous  systems  (DCCS)  (S child 
and  Lunze,  2008)  are  a  special  class  of  HDS  widely  used  in 
the  literature.  In  these  systems,  the  changes  in  discrete 
modes  are  achieved  by  discrete  control  commands,  e.g. 
opening  or  closing  a  switch. 

1.2  State  of  the  art 

Many  approaches  have  been  proposed  in  the  literature  for 
fault  diagnosis  of  DCCS.  They  are  generally  divided  into 
three  main  categories: 

•  approaches  for  the  diagnosis  of  parametric  faults, 

•  approaches  for  the  diagnosis  of  discrete  faults, 

•  approaches  for  the  diagnosis  of  both  parametric  and 
discrete  faults. 

In  parametric  fault  diagnosis  approaches,  (Cocquempot  et 
al .,  2004),  (Alavi  et  al.,  2011),  (Kamel  et  al.,  2012)  relations 
over  observable  variables  are  computed  in  order  to  generate 
residuals  sensitive  to  a  certain  subset  of  parametric  faults  in 
each  observable  discrete  mode. 

The  discrete  fault  diagnosis  approaches  are  divided  into 
three  main  groups.  In  the  first  group  (Rahiminejad  et  al., 
2012),  (Defoort  et  al.,  2011),  residuals  sensitive  to  the 
continuous  dynamics  in  each  discrete  mode  are  defined.  If 
unpredicted  change  occurs  due  to  the  occurrence  of 
unobservable  discrete  fault,  the  residuals,  defined  for  the 
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discrete  mode  before  the  fault  occurrence,  will  be  different 
of  zero  in  the  discrete  mode  after  the  fault  occurrence.  This 
change  of  residuals  values  from  zero  indicates  the 
occurrence  of  a  discrete  fault.  The  approaches  of  second 
group  (Bhowal  et  al .,  2007),  (Biswas  et  al. ,  2006),  describe 
in  each  normal  or  fault  discrete  mode,  continuous  dynamics 
as  the  rate  of  changes  of  continuous  variables.  These  rates 
are  considered  to  be  constant.  Transition  guards  are  defined 
as  linear  inequalities  based  on  continuous  variables  values. 
When  a  guard  is  satisfied,  its  corresponding  mode  transition 
is  enabled.  The  occurrence  of  a  fault  is  diagnosed  by 
determining  the  discrete  state  reached  due  to  specific  guard 
satisfaction.  In  the  methods  of  last  group  (Bayoudh  et  al ., 
2006),  a  set  of  residuals  is  defined  in  each  normal  or  fault 
discrete  mode.  Each  residual  is  characterized  by  three 
symbols:  0,  1  or  und  when  the  residual  value  is, 
respectively,  zero,  different  of  zero  and  undefined,  und 
represents  the  case  where  the  associated  residual  is  not 
defined  in  the  new  active  mode.  These  symbols  are  used  to 
distinguish  the  different  normal  and  fault  discrete  modes.  A 
discrete  fault  is  isolated  by  determining  the  current  discrete 
fault  mode  of  the  system. 

The  third  category  includes  few  approaches  for  the 
diagnosis  of  both  parametric  and  discrete  faults.  Some 
approaches  of  this  category  (Derbel  et  al.,  2009),  capture  the 
continuous  dynamics  by  integrating  the  occurrence  time  of 
events.  They  consider  that  the  occurrence  of  discrete  or 
parametric  faults  does  not  change  events  ordering  but  only 
alters  their  timing  characteristics.  Therefore,  a  discrete  or 
parametric  fault  is  diagnosed  when  predicted  events  occur 
too  late  or  too  early  or  they  do  not  occur  at  all  during  their 
predefined  time  intervals.  Other  methods  (Daigle  et  al., 
2010),  construct  temporal  causal  graphs  (TCG)  for  each 
normal  and  fault  discrete  mode  based  on  the  use  of  a  global 
hybrid  bond  graph.  When  measurement  deviations,  caused 
by  fault  occurrence,  are  observed  through  residuals,  TCG 
are  used  to  determine  the  effects  that  faults  will  have  on  the 
measurements  as  well  as  the  temporal  order  in  which  they 
deviate.  Then,  fault  signature  is  defined  for  each  fault  as  the 
qualitative  value  of  the  magnitude  and  the  first  non-zero 
derivative  change  which  can  be  observed  in  the  residuals.  In 
order  to  distinguish  parametric  from  discrete  faults,  the 
signatures  are  extended  by  adding  discrete  symbols 
indicating  abrupt  changes  from  zero  to  non-zero  or  from 
non-zero  to  zero.  In  (Louajri  et  al.,  2013),  an  approach 
based  on  a  diagnoser  with  hybrid  structure  is  developed.  It 
consists  of  three  parts:  the  discrete  diagnoser,  the  continuous 
diagnoser  and  the  coordinator.  The  discrete  diagnoser  is 
built  using  a  discrete  time  hybrid  automata  representing 
global  model.  It  exploits  the  information  extracted  from  the 
system  continuous  dynamics  to  get  rid  of  diagnosis 
ambiguity  due  to  the  system  behavior  abstraction.  The 
continuous  diagnoser  generates  residuals.  The  latter 
compare  the  measured  and  nominal  values  of  each 
continuous  variable  in  order  to  diagnose  the  parametric 


faults  in  each  discrete  mode.  The  information  about  the 
discrete  mode  is  provided  to  the  continuous  diagnoser 
thanks  to  the  information  extracted  from  the  discrete 
dynamics.  Finally,  the  coordinator  uses  the  decisions  issued 
from  the  discrete  and  continuous  diagnosers  in  order  to 
diagnose  faults  requiring  the  interaction  between  both 
diagnosers. 


1.3  Our  approach 

Fault  diagnosis  approaches  of  the  literature  do  not  scale  to 
HDS  with  a  large  number  of  discrete  modes  because  they 
achieve  fault  diagnosis  using  one  centralized  diagnosis 
module.  The  latter  is  built  using  a  global  model  of  the 
system.  Two  problems  are  arisen  -)  the  weak  robustness  in 
the  sense  that,  when  the  global  diagnosis  module  fails,  this 
may  bring  down  the  entire  diagnosis  task  and  -)  the  system 
global  model  can  be  too  huge  to  be  physically  constructed. 
Therefore  in  this  paper,  the  proposed  approach  of  (Fouajri  et 
al.,  2013)  is  developed  to  achieve  the  diagnosis  of 
parametric  and  discrete  faults  in  decentralized  manner  using 
several  local  hybrid  diagnosers.  The  latter  are  constructed 
without  the  use  of  a  global  model  of  the  system  but  only  the 
local  models  of  the  system  discrete  components  (Figure  1). 


The  paper  is  organized  as  follows.  In  section  2,  the  three  cell 
converter  system  is  described  and  modelled.  Section  3 
defines  the  steps  of  the  hybrid  diagnosis  construction.  In 
section  4,  a  simulation  for  the  three-cell  converter  is  used  to 
demonstrate  the  efficacy  of  the  approach.  A  conclusion  with 
the  future  work  ends  the  paper  in  section  5. 


Coordinator 

X 

Global  diagnosis  decision  DD 


Figure  1.  Decentralized  hybrid  diagnosis  structure  for  a 
HDS  composed  of  3  interacted  HCs. 


2.  Three  cell  converter  description  and 

MODELING 


2.1.  System  description 

In  order  to  illustrate  the  proposed  approach,  the 
decentralized  fault  diagnosis  of  three-cell  converters 
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(Shahbazi  et  al. ,  2013)  (Trigeassou,  2011),  depicted  in 
Fig.2,  is  achieved.  With  the  same  observability  used  in  the 
literature  (Defoort  et  al.,  2011),  (Uzunova  et  al.,  2012)  for 
the  three-cell  converter  diagnosis,  the  proposed  approach 
has  the  advantage  to  diagnose  (not  only  detect)  discrete  and 
parametric  faults  using  a  decentralized  structure. 

The  continuous  dynamics  of  the  system  are  described  by 
state  vector  X  =  [Vc±  Vc2  I]T ,  where  Vct  and  Vc2 
represent,  respectively,  the  floating  voltage  of  capacitors  C1 
and  C2  and  /  represents  the  current  flowing  from  source  E 
towards  load  ( R,L )  through  three  elementary  switching 
cells  Sj,j  E  {1,2,3}.  The  latter  represent  the  system  discrete 
dynamics.  Each  discrete  switch  Sj  has  two  discrete  states: 
Sj  opened  ( hq  =  0)  or  Sj  closed  ( hq  =  1),  where  hq  is  the 
state  discrete  output  of  Sj.  The  control  of  this  system  has 
two  main  tasks:  -)  balancing  the  voltages  between  the 
switches  and  -)  regulating  the  load  current  to  a  desired 
value.  To  accomplish  that,  the  controller  changes  the 
switches’  states  from  opened  to  closed  or  from  closed  to 
opened  by  applying  discrete  commands  ‘close’  or  ‘open’  to 
each  discrete  switch  Sj,j  E  {1,2,3}  (see  Fig.2).  Thus,  the 
considered  example  is  a  DCCS. 


Figure  2.  Three-cell  converter  discription  and 
decomposition 


2.2.  System  modeling  and  decomposition 

The  real  system  dynamic  evolution  of  three-cell  converter  is 
written  as  (Defoort  et  al.,  2011) 

(Vc  =  -hl-l  +  h2a-l 

1  q  Ci  q  C-L 

<Vc2  =  -h2q±I  +  h3q±I  (1) 

J  =  ~jl+h^\  Vc,  +  h2q  \  (Vc2  -VCl)+  h\\(E-  Vc2 ) 

As  shown  in  (1),  the  discrete  state  of  Sl9  represented  by  a 
real  discrete  output  /ij,  influences  the  dynamic  evolution  of 
Vcx  and  /.  The  discrete  state  of  S2 ,  represented  by  hq, 
impacts  the  dynamic  evolution  of  Vcx,  Vc2  and  /.  The 
discrete  state  of  S3,  represented  by  hq,  influences  the 
dynamic  evolution  of  Vcx  and  Vc2.  Thus,  the  three-cell 


converter  system  is  decomposed  into  three  interacted  HCs 
as  shown  in  Fig.2: 

•  HC1  is  composed  of  switch  S±  (Dcx),  Vc1  (Ccx)  and 
I  (Cc3). 

•  HC2  is  composed  of  switch  S2  ( Dc2 ),  Vcx  (Ccx), 
Vc2  (Cc2)  and  I  (Cc3). 

•  HC3  is  composed  of  switch  S3  ( Dc3 ),  Vc2  ( Cc2 )  and  I 
(Cc3). 

In  the  literature  (Defoort  et  al.,  2011),  (Uzunova  et  al., 
2012),  eight  faults  are  considered  for  the  diagnosis  of  the 
three-cell  converters  system  (Table  1). 


Table  1.  Faults  for  the  diagnosis  of  three-cell  converters 


Fault  types 

Fault  labels 

Fault  description 

Discrete  faults 

Fi 

S !  stuck  opened 

F2 

S1  stuck  closed 

Fi 

S2  stuck  opened 

F* 

S2  stuck  closed 

Fs 

S3  stuck  opened 

f6 

S3  stuck  closed 

Parametric 

faults 

Fv 

Change  in  the  nominal 
parameter  values  of  Ct  due 
to  Ci  ageing 

Fs 

change  in  the  nominal 
parameter  values  of  C2  due 
to  C2  ageing 

Labels  A^,  N2  and  N3  signify  the  normal  operating  modes 
for,  respectively,  HCl9  HC2  and  HC3. 

2.3.  Residuals  generation 

In  order  to  show  the  influence  of  each  discrete  component 
on  the  dynamic  evolution  of  each  continuous  component, 
(1)  is  rewritten  as  follows: 

Vc1  =  Vc\  +  Vc\ 

Vc2  =  Vcl  +  Vcl  (2) 

/  =  lc  +  71  +  /2  +  /3 

where  Vc\  =  -h\^l,  Vc\  =  h\^l,  Vc\  =  -h\^l, 
Vcl  =  h\±l,  ic  =  -~li,  i1  =  hl\vCl,  1 2  =  h2q-L(Vc2  - 
VCl),i3  =h3-L(E-Vc2). 

Vcl  represents  the  real  dynamic  evolution  of  Vcx  according 
to  the  discrete  state  of  ( Dct ).  Likewise,  Vcl ,  Vch 
i1,  i2 and  73  have  the  same  definition  as  Vc\.  ic  represents 
the  part  of  dynamic  evolution  of  /  which  does  not  depend  on 
the  discrete  state  of  any  switch. 

Similarly,  considering  that  the  parametric  faults  related  to 
the  load  ( R ,  L)  are  not  considered,  the  equations  system  for 
the  nominal  dynamic  evolution  of  system  components  can 
be  written  as: 
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f  vCl  =  Vc{  +  Vcl 

|  Vc2  =  Vcl  +  Vcl  (3) 

(j  =  ic  +  i1  +  i2  +  i3 

Vci  =  ^ ^ ^  = 

hq^rlJc  =tc  =  ~il  i1  =  h\-VCl,  I2  =  h2 -l  (Vc2  - 

Vcl,  /3  =  h3-L(E  —  Vc2)  where  hq,  hq  and  hq  are  the 
nominal  values  of  states  qlf  q2  and  q3  discrete  outputs  while 
C1  and  C2  are  the  nominal  values  of  C1  and  C2.  Based  on  (2) 
and  (3),  residuals  r1(  r2  and  r3  are  generated  as  follows: 


r2  =  (-h2q±+h2q±)I  +  (h3q±-h3q±)1 
r3  =  ^-hl)Vf+(h2-h2)^^ 
,  +m-h3q)^ 


(4) 


S^stuck-rtose,  fc^cj  denotes  the  set  of  fault  events 

(discrete  and  parametric)  that  can  occur  in  HC1.  The  set  of 
fault  events  contains  three  different  fault  types  or  modes 
indicated  by  the  fault  labels:  {FlfF2,F7}.  The  set  of  labels 
for  HC1  is  SP1  =  {  Nlf  Flf  F2>  F7}  . 

51:  Q1  X  X1  ->  Q1:  is  the  state  transition  function.  A 
transition  51(q1,  e)  =  q1+  corresponds  to  a  change  from 
state  q1  to  state  q1+  after  the  occurrence  of  event  e  E  X1; 

X 1  =  { VclfVc2,I }  is  a  finite  set  of  continuous  variables 
associated  to  S1 ; 

flux1-.!}1  x  X1  ->  Mn={X1,X1}:  is  a  function 

characterizing  temporal  evolution  X1  and  nominal  evolution 
X1  of  continuous  variables  X1  in  each  discrete  state  ql, 
where  i1  =  [vd  Vc\  l1]* ,  X1  =  [Vcl  Vcl  T]r; 

Init 1  c  Q1  x  X1  =  S±0  (c/o):  is  the  set  of  initial 
conditions. 


In  order  to  show  the  influence  of  each  discrete  component 
on  the  residuals,  (4)  is  rewritten  as  follows: 

ri  =  ri  +  ri 

r2  =  r|  +  r|  (5) 

r3  =  r3c  +  rl  +  r32  +  r| 

where  r/  =  (— hq  ^-)  /  =  {vc\  —  Vc^,  r?=(vcl  — 

Vcl'),  r£=(Vcl  -  Vcf),  r2  =  (Vc%  -  lYf),  rl  = 
(j1  —  i1  r32  =  {i2  —  i 2),  r3  =  (/3  —  /3)  and  r3c  = 

(4  -  4)  =  0. 

2.4.  Hybrid  automata  construction 

Hybrid  automata  A1  characterizing  the  hybrid  dynamics  of 
HC1  is  defined  by  the  tuple  (see  Fig.4  and  Fig.4): 

A1  =  (Q\ X1,  SP \  S1,  X1,  flux \  r1,  Init1)  (6) 

where, 

Q1={S10  ( S1  opened),  (5!  closed),  5X50  (S1  stuck 

opened),  5X5C  (5X  stuck  closed)}:  is  a  finite  set  of  discrete 
states  (discrete  modes)  of  S1.  The  output  of  state  ql  is 
characterized  by  real  discrete  output  vector  hq  E 
(0  (when  S1  is  opened),  1  (when  S1  is  closed)}  and  nominal 
discrete  output  vector 

hq  =  {0  (when  S±  have  to  be  opened),  l(when  S±  have  to  be 
closed) } .  At  normal  discrete  mode  (state)  h\  =  hq  while  in 
faulty  mode  hq  =£  /ij; 

X1  =  X\  U  X1  is  the  event  set  of  S1 .  It  includes  observable 
events  corresponding  to  control  command  events  X1  = 
(close  5X),  OS^open  5^)}  and  unobservable  events  X1 
including  fault  events,  xj  =  {5  x_stuck_oipen. 


Name  of  state 


h'q,hq 


Vcl 

Vcl' 

rl  =  Vc\  -Vcl 

Vcl 

x\  = 

Vcl 

II 

I 

II 

I1 

i1 

II 

L  -1 

i  SP/  i 

Figure  3.  Hybrid  state  of  A1  for  HC1. 


Figure  4.  Hybrid  automata  A1  for  HC1. 


r1  =  {rf,  r2  ,  731}:  is  a  set  of  residuals  associated  to  HC 
Since  V c2  does  not  belong  to  HC ±,  therefore,  Vc2  =  0, 
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V  c\  =  0.  Thus  r2  is  equal  to  zero. 

Hybrid  automata  A2  and  A3  for  HC2  and  HC3  are 
constructed  by  the  same  manner. 

2.5.  Motivation  to  use  the  considered  residuals 

Let  us  consider  the  occurrence  of  a  fault  of  type  Fs,  e.g. 
S3  stuck  opened.  When  the  controller  sends  control 
command  'CS3'( close  S3),  S3  remains  in  its  stuck-opened 
mode  ( [hq  =  1  and  h3q  =  0).  The  occurrence  of  a  fault  of 
type  Fs  impacts  at  the  same  time  r2  and  r3  (h3  =  1  and 
hq  =  0)  while  it  does  not  impact  ((hq  =  hq),  (hq  =  hq) 
and  ( C1  =  C±)),  see  (4).  Therefore,  there  is  no  delay  of  the 
influence  of  the  fault  occurrence  on  the  sensitive  residuals, 
e.g.  r2  and  r3.  Moreover,  there  is  no  fault  propagation  from 
one  residual  to  another  one,  from  r2  or  r3  towards  rt . 

3.  Three  cell  converter  diagnosis 

3.1.  Global  fault  signature  construction 

A  qualitative  signature  is  constructed  by  generating 
continuous  and  discrete  symbols  from  residual  values. 
Continuous  symbols  CSfr^  E  (0,  — ,  +}  represent  the 
qualitative  abstraction  of  residual  values  into 
stable/increasing/decreasing  ones : 

•  rf :  rt  (t)  belongs  to  the  nominal  interval; 

•  n  :  ri  (0  *s  below  the  nominal  interval; 

•  ri+ :  ri  (0  is  above  the  nominal  interval. 

The  occurrence  of  a  discrete  fault  exhibits  an  abrupt  change 
in  the  continuous  dynamics  due  to  unpredicted  change  in 
Dcj  discrete  mode.  This  change  is  characterized  by  the 

absence  (  hJq  =  0  while  hq  =  1)  or  the  addition  ( hq  =  1 
while  hq  =  0)  of  associated  term  e.g.,  — .  On  the  other 

hand,  parametric  faults  due  to  the  ageing  effect  cannot  cause 
this  abrupt  change  with  a  finite  change  in  magnitude.  In 
fact,  they  are  indicated  by  a  progressive  abnormal  change  of 
the  parameter  value.  In  order  to  take  into  account  this 
discriminative  information,  discrete  symbols  DSfji)  are 
added  for  the  abstraction  of  each  residual  rt  in  order  to 
distinguish  between  parametric  and  discrete  faults  as 
follows: 

•  PC{  =  +Val:  denotes  an  abrupt  positive  change  in 
residual  ri  due  to  a  discrete  fault  caused  by  Dcj. 
+Val  is  equal  to  the  absolute  value  of  the  term 
associated  to  hq', 

•  NC(  =  —Val:  denotes  an  abrupt  negative  change  in 
residual  rt  due  to  a  discrete  fault  caused  by  Dcj ; 

•  UCf.  denotes  that  there  is  no  observed  abrupt  change 
in  residual  rt . 


•  A  fault  signature  Sigq  at  global  discrete  state  q  is  the 
combination  of  continuous  and  discrete  symbols  of 
the  different  residuals  as  follows: 

Sigq  =  (r“(ri  \  DSdn  ))&■■■  &(rnS(rn  \  DS(rn  ))  (7) 

3.2.  Local  fault  signature  construction 

Each  discrete  state  qJk  of  Aj  generates  a  fault  signature  sigJk 
as  a  guard  over  residuals  r7  calculated  in  this  discrete  state 
as  follows: 

sigl  =  (V  CS{r>1  \  DS(r -/)j  & ...  &  (rj  cs(r"  \  D5(rn;'))  (8) 

Based  on  (5),  we  can  write: 

n=*i~  *i=  -  x})+...+  (x[  -  x[)  =  r/+. . .  +r/ 

If  [(x/  —  x[)  =  r(J]  ^  0,  it  means  that  the  other  parts  of 
residual  rt  are  equal  to  zero  (one  fault  can  be  occurred  at  the 
same  time).  In  this  case,  rt  =  r-; .  Hence,  rt  will  have  the 
continuous  and  discrete  symbols  of  r-  .  Thus  (8)  is  rewritten 
as  follows: 

si9k  =  (rP(ri\DS(r1)^8i...&{r^rn),DS(rn))  (9) 

By  comparing  (8)  and  (9),  we  can  notice  that  sigk  becomes 
equivalent  to  the  global  fault  signature  Sigq. 

3.3.  Local  hybrid  diagnoser 

The  objective  of  local  hybrid  diagnoser  Dj  is  to  detect  and 
isolate  the  occurrence  of  parametric  and  discrete  faults 
affecting  the  dynamics  of  hybrid  component  HCj.  Dj  is  built 
based  on  the  local  model,  A7,  of  HCj.  Each  state  of  £>;, 
denoted  zk ,  is  of  the  form  shown  in  Fig. 5. 


_ 4 _ 

Model  states:  Ql_ 

_ xj_ _ 

SPJ 

Figure  5.  State  of  local  hybrid  diagnoser  Dj  of  HCj. 

Local  hybrid  diagnoser  D1  of  HC1  is  depicted  in  Fig. 6.  It  is 
constructed  from  hybrid  automata  A1  of  Fig.4. 

is  constructed  as  follows: 

•  Initial  state  characterized  by  ( Ql ,  X1,  SP1),  is 
composed  of  the  following  A1  states:  ql  (A1  initial 
state),  q\  reached  from  q\  by  the  occurrence  of  a  fault 
event  ‘S^stuck-Open’  (fault  of  type  /q)  and  ql  reached 
from  ql  due  to  the  occurrence  of  a  fault  event  ‘fc^c^ 
(fault  of  type  F7).  Thus,  Ql  is  equal  to  {ql,  q3,  q)}.  SP 1 
gathers  the  normal  and  fault  labels  associated  to  the  states 
belonging  to  Ql.  Therefore,  SP 1  is  equal  to  {N±,  F1,  f7}. 


269 


Annual  Conference  of  the  Prognostics  and  Health  Management  Society  2014 


Finally,  X1  gathers  of  all  the  states  ql  of  Q\.  Since 
states  q\  and  q7  are  reached  from  q\  due  to  the 
occurrence  of  unobservable  event  (a  fault),  X\,  X3  and 
X7  are  equivalent  and  equal  to  [0  0  0]r  (see  Fig.4). 

•  The  states  of  Dx  reached  due  to  the  occurrence  of  each 
control  command  event  observed  by  HC1  are  computed. 
Since  D±  initial  state  is  5X0,  control  command  0S1  will 
not  change  Dx  state  z\.  The  event  CS1  transits  D±  from 
z\  to  z\,  characterized  by  X1,  SP 1).  Q\  is  equal  to 
all  the  states  reached  from  Q\  due  to  the  occurrence  of 
CS1.  Thus,  Q\  is  equal  to  {q\,q\,  ql }  (see  Fig.4). 
Moreover,  all  the  states  of  A1  reached  from  Q\  due  to  the 
occurrence  of  unobservable  event  are  added  to  Q\. 
Therefore,  Q\  is  equal  to  {q\,q\,  ql,  ql).  SP 1  is  equal 

to  { Nl9  F1?  F7,  F2 }  while  X1  is  equal  to  0  (see 

Fig.4). 

•  Fault  signatures  are  generated  for  each  D\  state  thanks  to 
the  continuous  dynamic  evolution  in  each  discrete  state 
of  Ql.  In  the  initial  Dt  state,  z\,  the  continuous  dynamic 
evolution  in  any  state  of  Q\  does  not  evolve.  Therefore, 
their  associated  residuals  are  equal  to  zero  leading  to 
obtain  the  fault  signature  Sigh  (see  Table  2).  In  z\,  the 
continuous  dynamic  evolution  of  the  states  belonging  to 
Ql  will  allow  to  generate  four  fault  signatures  as  we  can 
see  in  Fig. 6.  They  allow  to  detect  and  isolate  discrete  and 
parametric  faults  F1  and  F7  as  follows.  q\  of  A1  (reached 
due  to  the  occurrence  of  fault  of  type  Ft)  generates  local 
fault  signature  Sig\. 

44 

Sigh  =  (ri1~’^)&(r210'UC2)&(r31+’T‘)  (see  the  values 
of  local  residuals  in  q\  of  Fig.4).  As  explained  in  subsection 
3.2,  local  fault  signature  Sig\  is  equal  to  global  fault 

44 

signature. 

Sigq  =  (rf^.DSfaj)  &(r2cs<r2,,DS(r2))& 

(r3cs™,DS(r3))  =  (rr,^)&(r2°,UC2)&(r+,^). 

This  global  signature  is  used  as  transition  to  isolate  the 
occurrence  of  a  fault  of  type  F1.  Same  reasoning  can  be 
applied  for  the  other  fault  signatures.  To  overcome  the  noise 

problem,  the  values  of  comparison  (e.g.,— )  are  replaced  by 

C2 

the  intervals  corresponding  to  the  selected  confidence  level. 
These  intervals  are  calculated  using  Z-test  in  order  to 
determine  the  thresholds  of  each  value. 

Same  reasoning  can  be  followed  for  the  construction  of  the 
other  states  of  D1 . 

It  is  worth  pointing  out  that 

Sigq  =  Sig1^  =  (rf , ^ &(r2°, UC2)& (r3+ , means  that 
the  three  conditions  have  to  be  satisfied  in  order  to  enable 


the  corresponding  transition. 


Figure  6.  Local  hybrid  diagnoser  D1  of  HC1. 

Table  II  shows  the  local  fault  signatures  (equivalent  to  the 
global  fault  signatures)  used  by  to  achieve  its  local 
diagnosis. 


Table  2.  Local  fault  signatures  generated  due  to  the 
_ occurrence  of  faults  in  HC^. _ 


SP1 

Local 

signature  name 

Equivalent  global  fault  signatures 

Pi 

Sig  i 

{r'%)Hr?’ucMr3+'ir) 

p2 

Sigl 

(r1+,^)&(r20,(/C2)&(r3-,^i) 

Py 

Sigh 

(rf,  UC,  )&(r2°,  UC2)&(rl  UC3) 

Si-972 

(r,+,UC,mrlUC2)HrlUC3) 

Nt 

Sigh 

(r1°,C/C1)&(r2°,yC2)&(r‘,,yC3) 

The  other  diagnosers  D2  and  D3  for  FIC2  and  HC3  can  be 
constructed  similarly  as  for  Dx .  D2  is  sensitive  to  discrete 
faults  F3  and  F4  and  to  parametric  faults  F7  and  F8,  while  D3 
is  sensitive  to  discrete  faults  F5  and  F6  and  to  parametric 
fault  F8.  The  occurrence  of  parametric  fault  F7  (respectively 
F8)  is  detected  intrinsically  by  Dx  and  D2  (respectively  D2 
and  D3). 

3.4.  Coordinator  construction 

The  system  decomposition  achieved  by  the  proposed 
approach  allows  each  local  hybrid  diagnoser  to  diagnose 
faults  that  can  occur  in  its  corresponding  hybrid  component. 
In  order  to  obtain  a  decentralized  diagnosis  performance 
equivalent  to  a  centralized  diagnoser,  a  decision  coordinator 
is  defined.  It  generates  a  global  diagnosis  decision  by 
merging  local  diagnosis  decisions  provided  by  local  hybrid 
diagnosers.  Let  us  denoted  F1,  F2  and  F3  the  faults  that  can 
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occur,  respectively,  in  HC1,  HC2  and  HC3.  F 1  E  [Flf  F2,  F7}, 
F2  E  {F3,F4,  F7,  F8}  and  F3  E  (F5,F6,  F8},.  Global  diagnosis 
decision  DD  is  computed  as  follows: 


In  order  to  highlight  the  efficiency  of  the  diagnoser,  the 
simulations  take  into  account  the  set  of  faults  defined  in 
Table  1  for  the  three-cell  converter. 


•  Dr  diagnoses  with  certainty  the  occurrence  of  a  fault 
of  type  F1  through  the  global  fault  signature  Sigq. 
D2  cannot  diagnose  with  certainty  the  occurrence  of 
this  fault  because  it  does  not  belong  to  its  associated 
HC2.  D3  cannot  diagnose  with  certainty  the 
occurrence  of  this  fault  because  it  does  not  belong  to 
its  associated  HC3.  Therefore,  the  global  diagnosis 
decision  will  be  DD  =  F1. 

•  Global  fault  signature  Sigq  corresponds  to  a  fault  of 
type  F1  or  of  type  F 2  ( F7 ).  Thus,  global  diagnoser 
DD  will  be  F1  or  F2 .  Both  D±  and  D2  are  sensitive  to 
this  fault  signature,  therefore  Dx  declares  F7;  and  D2 
declares  F7.  In  order  to  obtain  a  decentralized 
diagnosis  decision  equivalent  to  the  global  one, 
global  diagnosis  decision  DD  will  be  equal  to 
(F1  or  F2)  =  F7. 

•  Table  3  shows  global  diagnosis  decision  DD.  A  local 
diagnoser  declares  ‘nothing’  when  it  cannot  confirm 
the  occurrence  or  the  non-occurrence  of  a  fault. 

Table  3.  Global  diagnosis  decision  DD  for  Three  Cell 


Converter  . 


cases 

Local 

Local 

Local 

Global 

decision 

DD 

diagnoser  D1 

diagnoser  Dz 

diagnoser  D3 

1 

W, 

N* 

N* 

N 

2 

F1 

N2  or  Nothing 

N2  or  Nothing 

F1 

3 

N-l  or  Nothing 

F2 

N3  or  Nothing 

F2 

4 

F 1 

F2 

N3  or  Nothing 

F1or  F2 

5 

N-l  or  Nothing 

F2 

F3 

F2ot  F3 

6 

N-l  or  Nothing 

N2  or  Nothing 

F3 

F3 

1 

Nothing 

Nothing 

Nothing 

Nothing 

3.5.  Identification  of  parametric  faults 

When  one  of  parametric  faults  is  diagnosed,  its  real  value 
needs  to  be  identified.  As  an  example,  for  parametric  fault 
of  type  F7  related  to  Cl9  the  real  value  of  the  latter  is 
identified  based  on  its  corresponding  residual  as  follows: 

r>  =  (-s+s)=c'=Ai  <10» 

The  same  reasoning  is  applied  to  identify  the  real  value  of 
capacitor  C2  in  case  of  fault  of  type  F8  related  to  C2 . 

4.  Experimentation  and  obtained  results 

In  order  to  evaluate  the  proposed  approach,  simulations 
were  carried  out  for  the  three-cell  converter  using  Matlab- 
Simulink™  environment  and  Stateflow™  toolbox.  The 
parameters  used  in  these  simulations  are: 

F  =  60V,  C1  =  C2=  40 gF,R  =  200/2,  L  =  0.1  H. 


Discrete  controller  commands  are  assured  by  a  pulse  width 
modulation  (PWM)  signal  (Defoort  et  al .,  2011).  Fig.7 
depicts  the  control  of  three  switches  S±  ,S2  and  S3 .  When  the 
triangular  signal  is  below  the  reference  signal  (ref  in  Fig.7), 
the  associated  switch  is  controlled  to  be  opened.  When  the 
triangular  signal  is  above  the  reference  signal,  the  associated 
switch  is  controlled  to  be  closed.  This  sequence  of  control  is 
periodic  with  a  period  of  TPWM  =  0.025. 


S, 


S2 


S>- 

°0  0.005  0.015  0.025  0.035  t(s) 

Figure  7.  PWM  for  control  of  three  switches  51?  S2and  S3. 

4.1.  Normal  conditions  scenario 


Fig. 8  depicts,  respectively,  the  signals  of  floating  voltages 
Vct  and  Vc2  and  the  current  /.  These  signals  correspond  to 
the  normal  conditions.  Moreover,  one  can  see  in  Fig. 8  that 
Vcx  (respectively  Vc2)  has  a  periodic  signal  corresponding 

to  load  and  unload  of  capacitor  Ct  (respectively  C2)  around 

£ 

the  mean  value  Vclref  =  -=  20V  (respectively Vc2ref  = 

2  E 

—  =  4 OF)  and  that  the  current  /  remains  constant  in  the 
region  of  its  reference  value  (0.15A). 


Fig. 9  shows  the  real  and  nominal  dynamic  evolution 

ofVc^Vc^ndVcJ,  Vc2(  Fc2and  Vc2)  and  /(/and/).  We 
can  notice  that  the  curves  representing  the  real  and  nominal 
dynamic  evolutions  are  superposed.  Consequently, 
residuals  rl9  r2  and  r3  are  equal  to  zero  in  these  conditions. 


0  0,1  0,2  0,3  0,4  0,5  0,6  0,7  0,8  0,9  1  t(s) 


Figure  8.  Real  signals  corresponding  to  Vcl9  Vc2  and  /  in 
normal  conditions  . 
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Figure  9.  Real  and  nominal  dynamic  evolution  of  Vcl9  Vc2  and  /  in  normal  conditions. 


4.2.  Faulty  conditions  scenario 

The  test  scenario  is  generated  as  follows  (see  Fig.  10).  Each 
fault/,  belonging  to  one  of  the  fault  labels  of  Table  1,  is 
generated  starting  at  time  tsf  and  ending  at  time  tef.  Then,  the 
system  returns  to  normal  operating  conditions  before 
generating  a  new  fault  for  a  certain  time.  Parametric  faults 
of  types  F7  and  F8  are  simulated  by  changing  gradualy  the 
real  values  of  Cl9  respectively  C2,  in  positive  or  negative 
direction  using  a  ramp  signal.  V cx  ,  Vc2  and  /  simulated 
signals  including  these  faults  are  represented  in  Fig.  11. 

One  can  see  in  Fig.  11  that  Vcx  (respectively  Vc2)  has  lost 
the  periodic  aspects  in  the  case  of  fault  and  that  the  current  / 
has  become  nonconstant  in  the  region  of  its  reference  value. 


r1?  r2,  r3  are  represented  in  Fig.  12  and  Fig.  13.  As  expected, 
rx  is  sensitive  to  the  faults  of  types  Fl9  F2,  F3,  F4  and  F7,  r2is 


Figure  10.  Time  of  apperance,  injection,  of  faults  during 
the  simulation  of  three  cell  converter. 
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Figure  11.  Real  signals  of  V cl9  Vc2  and  /  in  faulty  and 
normal  conditions. 


sensitive  to  the  faults  of  types  F3,  F4,  F5,  F6  and  F8  while  r3 
is  sensitive  to  the  faults  of  types  Fl9  F2,  F3,  F4,  F5  and  F6. 

Fig.  14,  Fig.  15,  Fig.  16  and  Fig.  17  show,  respectively,  local 
decision  (5FX)  of  diagnoser  Dl9  local  decision  ( SP2 )  of 
diagnoser  D2 ,  local  decision  (SP3)  of  diagnoser  D3  and 
global  decision  (SP). 

The  first  local  diagnoser  Dx  is  sensitive  to  faults  of  types  Fl9 
F2  and  F7  (diagnosis  with  certainty  their  ocurence),  the 
second  local  diagnoser  D2  is  sensitive  to  faults  of  types  F3, 
F4,  F7  and  F8  while  the  third  local  diagnoser  D3  is  sensitive 
to  faults  of  types  F5,  F6  and  F8.  We  can  conclude  that  the 
global  decision  indicates  with  certainty  the  occurrence  of 
each  of  the  generated  faults.  The  diagnosis  delay 
corresponds  to  the  time  when  the  system  is  in  a  discrete 
fault  is  due  to  residues  that  are  silent  in  some  discrete  state. 
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Figure  12.  Residuals  corresponding  to  generated  discrete 
faults  of  Fig.  10. 
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Figure  13.  Residuals  corresponding  to  generated 
parametric  faults  of  Fig.  10. 
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Figure  17.  Global  diagnosis  decision  issued  by  the  coordinator. 


4.3.  Normal  conditions  with  noises  in  parameters 
scenario 

Diagnosis  algorithms  should  be  tested  and  evaluated  on  real 
systems  with  practical  significance.  In  these  systems,  factors 
such  as  noise  make  diagnosis  challenging.  Therefore,  there 
is  a  need  to  evaluate  the  robustness  of  the  diagnosis 
algorithms  for  different  fault  and  noise  magnitudes. 


Accurate  simulation  models  of  the  system  are  required  for 
this  purpose.  Further,  it  is  important  to  execute  the  diagnosis 
algorithms  on  systems,  where  model  uncertainty  is  always 
present,  and  complicates  the  diagnosis  task.  In  order  to 
examine  the  robustness  of  our  approach,  a  parametric  noise 
(see  for  example  Fig.  18),  applied  on  parameters,  is  used. 
From  an  electrical  point  of  view,  the  resistors  are  the  most 
disturbing  element  in  tree  cell  converter  systems.  For  this 
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reason,  we  simulated  noise  on  signal  resistance. 


Figure  18.  Noise  added  to  resistance  R  in  the  converter. 

In  order  to  take  into  account  the  noises  in  R ,  the  residuals  of 
(4)  is  written  as  follows: 


=  (-%  i +h«rJ 1 + i 1 

<  r2  =  (~h2q  +K  A)  I  +  (h^2  -hi  A)  / 

r3  =  (-r  +  Rb)L  +  (h*  _/li)^  + 

Where  R  is  the  nominal  value  of  R  without  noises  while  Rb 
is  the  real  value  of  R .  The  latter  corresponds  to  the  nominal 
value  of  R  with  noises. 

rl9  r2,  r3  are  represented  in  Fig.  19.  As  expected,  rx  and  r2 
are  not  sensitive  to  this  perturbation  in  normal  conditions  ( R 
does  not  influence  the  dynamic  evolution  of  V cx  and  Vc2. 
While  r3  is  impacted  by  this  noise.  It  changes  between 
— OAA/s  and0.4A/s. 
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Figure  19.  Set  of  residuals  with  noise  corresponding  to  the 
normal  conditions. 


Ideally,  any  non-zero  residual  value  implies  a  fault,  which 
should  trigger  the  fault  isolation  system.  Therefore, 
statistical  techniques  are  required  for  reliable  fault  detection. 
The  fault  detection  system  is  based  on  a  Z-test  that  uses  the 
estimated  variance  of  the  residuals  and  a  pre-specified 
confidence  level  to  establish  the  significance  of  observed 
nonzero  residuals.  To  cope  with  noise,  we  compute  the 
mean  and  the  variance  at  different  time  points  (Biswas  et 


al. ,  2003).  The  Z-test  is  a  statistical  inference  test  employed 
to  establish  the  signification  of  the  deviation.  It  requires  the 
mean  and  standard  deviation  of  the  population,  and  the 
mean  and  size  of  the  samples.  These  values  are  estimated 
using  sliding  windows  over  the  residual  for  a  variable.  A 
small  sliding  window  of  size  W1  =  5  samples,  is  used  to 
estimate  the  current  mean  fir.  (  t)  of  the  residual  rt  related  to 
the  variable  xt : 

Mn(0  =  ^2r=t-w1+in(0  (13) 

We  suppose  the  mean  of  the  population  is  equal  to  zero, 
since  the  residual  should  be  zero  when  the  system  is  free  of 
faults.  We  compute  the  variance  from  data  history  of  the 
nominal  residual  signal  over  a  window  W2  proceeding  W1 
as  an  estimate  of  the  true  variance: 

=  —J^v=t-w2-W1  +  iri(V^  (14) 

°W(t)  =  ^-tv™-W2-Wl+ 1  (rA)  -  Mr, A  ))2  (15) 

The  size  of  W2  must  contain  enough  of  measurements  in 
order  to  estimate  correctly  the  residuals  ’  mean  and  variance 
in  the  normal  operating  conditions  and  therefore  to  reduce 
the  rate  of  false  alarms.  The  size  of  must  also  be  selected 
as  a  tradeoff  between  the  delay  of  fault  detection  and  the 
rate  of  false  alarms.  The  size  of  VIZ2 ,  respectively  W1,  is 
chosen  experimentally  to  be  equal  to  25,  respectively  5, 
measurements. 


Since  the  distribution  of  residuals  mean  is  supposed  to 
follow  the  normal  distribution,  a  confidence  level,  a,  is 
defined  by  determining  the  bound  Mr;  ]  within  which 
fir.  ( t)  is  considered  to  correspond  to  normal  operating 
conditions.  \^.r  Mr;  ]  is  defined  using  Z-test  table  and  the 
approximation  ar  : 


  zvi°ri 


(16) 


Mr, 


(17) 


For  a  equal  to  0.95,  zv.  and  are  equal  to,  respectively, 
-1.64  and  1.64. 


The  Z-test  is  employed  in  the  following  manner: 

\irt  <  Mr;  <  Mr;  =>  No  fault 
Otherwise  =>  Fault 

Fig. 20  depicts  mean  of  residuals  Mr3  and  the  negative  and 
positive  thereshold  of  this  residual.  The  mean  and  true 
variance  of  residual  rx  and  r2  are  equal  to  zero.  Thus  its 
thershold  is  also  equal  to  zero  (Mri,  Mr^  and  Mrx 
respectively,  fir. ,  Mr2and  Mr2are  superposed). 
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Figure  20.  Set  of  residuals  and  thersholds  with  noise 
corresponding  to  the  normal  conditions. 

In  case  of  fault,  Table  4  is  used  to  achieve  a  local  diagnosis 
of  D±. 


Table  4.  Local  fault  signatures  generated  due  to  the 
occurrence  of  faults  in  HCj  in  case  of  parametric  noise. 


SP1 

Local 

signature 

name 

Equivalent  global  fault  signatures 

Pi 

Sig\ 

()r  +  <  <  ^+^i)  &0*r2  <  Pr2  <  m£.)& 

Gr+^  <5r+A) 

Pi 

sigl 

+  Pn  <Pr2<^  +Pn)  &(Pr2  <  Pr2  < 

{ZXT  +  Pr3<Pr3<ZXT+Pr2) 

Pi 

Si-971 

(Mr,  >  MrJ&fe  <  Mr?  <  Mr?)&(MT,  <  Mr,  <  Mr,) 

Si-972 

(Mr,  <  Mr,)&(Mr,  <  Mr,  <  Mr,)&(Mr,  <  Mr,  <  Mr,) 

Ni 

Sigh 

(Mr,  <  Mr,  <  M?J&(Mr2  <  Mr2  <  M?2)&(Mr3  <  Mr3 

<o 

The  other  diagnosers  D2  and  D3  for  HC2  and  HC3  can  be 
constructed  similarly  as  for  . 


4.4.  Faulty  conditions  with  parameters  perturbation 

In  order  to  evaluate  the  proposed  approach  in  case  of  noise, 
another  scenario  of  fault  is  generated  (see  Fig.21).  The 
corresponding  iir± ,  (iT2 ,  /2r3  for  this  senario  are  represented  in 
Fig.22.  and  Fig.23.  In  this  case,  noises  are  observed  only  in 
r3  at  normal  and  faulty  conditions  (see  zoom  in  Fig. 24).  As 
we  said  before,  only  r3  is  impacted  by  noises  since  the  noisy 
parameter  R  is  included  only  in  dynamic  evolution  /  of  I  (see 
(1)).  To  overcome  this  problem,  a  threshold  is  defined  for 
each  residual  using  Z-test.  These  thresholds  are  used  during 
the  fault  detection  and  isolation  in  order  to  avoid  the  false 
alarms  as  well  as  the  fault  missed  detection  caused  by 
noises. 
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Figure  22.  Residuals  corresponding  to  generated  discrete 
faults  of  Fig.21  in  case  of  noise. 
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Figure  23.  Residuals  corresponding  to  generated 
parametric  faults  of  Fig.21  in  case  of  noise. 


Figure  24.  Zoom  of  residuals  signals  with  noise 
corresponding  to  normal  and  faulty  conditions. 


Figure  21.  Time  of  apperance,  injection,  of  faults  during 
the  simulation  of  three  cell  converters  with  noise. 


Fig. 24,  Fig. 25,  Fig. 26  and  Fig. 27  show,  respectively,  local 
decision  (57\)  of  diagnoser  Dl9  local  decision  (SP2)  of 
diagnoser  D2 ,  local  decision  (SP3)  of  diagnoser  D3  and 
global  decision  ( SP ).  The  first  local  diagnoser  is 
sensitive  to  faults  of  types  Fl9  F2  and  F7  (diagnosis  with 
certainty  their  ocurence),  the  second  local  diagnoser  D2  is 
sensitive  to  faults  of  types  F3,  F4,  F7  and  F8  while  the  third 
local  diagnoser  D3  is  sensitive  to  faults  of  types  F5,  F6  and 
F8.  We  can  conclude  that  the  global  decision  indicates  with 
certainty  the  occurrence  of  each  of  the  generated  faults 
regardless  of  the  existence  of  noise. 
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Figure  25.  Local  decisions  ( SP \)  of  D1  in  case  of  noise. 


Figure  26.  Local  decisions  ( SP2 )  of  D2  in  case  of  noise. 


Figure  27.  Local  decisions  ( SP3 )  of  D3  in  case  of  noise. 


0.1  0.15  0.2  0.25  0.3  0.35  0.4  0.45  0.5  0.55  0.6  0.65  0.7 


0.8 


0.9  0.95  1  t(s) 


Figure  28.  Global  diagnosis  decision  issued  by  the  coordinator  in  case  of  noise. 


5.  Conclusion 

In  this  paper,  a  decentralized  hybrid  diagnosis  approach  for 
discretely  controlled  continuous  systems  is  proposed.  The 
elaboration  of  this  approach  is  motivated  by  the  capacity  of 
the  hybrid  models  to  represent  intrinsically  the  interactions 
between  the  continuous  and  the  discrete  dynamics  of  a 
system. 

The  originality  of  this  work  is  the  exploitation  of  the  system 
modularity  in  order  to  reduce  its  complexity  as  well  as  the 


explosion  in  the  number  of  its  discrete  states.  To  achieve 
that,  the  diagnosis  task  is  accomplished  by  a  set  of  local 
hybrid  diagnosers.  Each  of  the  latter  is  responsible  of  the 
diagnosis  of  a  specific  part  of  the  system.  These  local  hybrid 
diagnosers  are  built  without  the  use  of  the  system  global 
model  but  only  local  models.  The  decisions  of  the  local 
hybrid  diagnosers  are  merged  using  a  coordinator  in  order  to 
obtain  a  diagnosis  performance  equivalent  to  the  one  of  a 
centralized  diagnosis  structure. 

In  the  future  work,  this  approach  will  be  applied  to  a  real 
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three-cell  converter.  Then,  it  will  be  developed  to  consider 

multiple  and  adjacent  faults  in  a  more  general  class  of 

hybrid  dynamic  systems. 
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